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The recent reports on high-temperature superconductivity above 190 K in hydrogen sulfide at 
200 GPa pressure, exceeding all previously discovered superconductors, has greatly invigorated in¬ 
terest in dense hydrogen-rich solids. Here, we investigate a possible way to optimize the critical 
temperature in these compounds by using first-principles linear-response calculations. We construct 
hypothetical alchemical atoms to smoothly interpolate between elements of the chalcogen group and 
study their bonding and superconducting properties. Our results show that the already remarkable 
critical temperatures of H 3 S could be improved even further by increasing the ionic character of the 
relevant bonds, i.e., partially replacing sulfur with more electronegative elements. 


The recent experimental report of superconductivity 
with a critical temperature (T c ) exceeding 190 K in dense 
hydrogen sulfide [I] could be one of the most important 
discoveries in the field of superconductivity, if confirmed. 
Although x-ray scattering data are lacking due to the 
extreme pressures involved (p ~ 200 GPa), there is a 
widespread consensus that the superconducting samples 
are composed of a sulfur-trihydride (H 3 S) bee phase pre¬ 
dicted by Duan et ai [2]. Theoretical studies [3H9] and 
isotope effect measurements uni unambiguously indicate 
that, unlike the high-T c iron pnictides and cuprates, H 3 S 
is a conventional BCS m superconductor. 

Phonon-mediated superconductors are accurately de¬ 
scribed by the Migdal-Eliashberg theory Hu eng; their 
T c is well approximated by the Alien-Dynes formula [14] , 


rji ^ln 
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where A is the electron-phonon ( e-p ) coupling constant, 
co>i n is the logarithmic average of phonon frequencies, 
and n* is the Coulomb pseudopotential. Thus, Migdal- 
Eliashberg theory indicates a clear path to optimize T c , 
i.e., looking for materials with high characteristic phonon 
frequencies (cji n ) and strong e-p coupling constants A. 

In a seminal paper that is now almost 50 years old, 
Ashcroft noticed that these conditions could be realized 
in ultra-dense metallic hydrogen HE yet only the recent 
progress in experimental techniques for high-pressure re¬ 
search and improved methods for crystal-structure pre¬ 
diction have made this hypothesis viable [16 . In the last 
few years, there has been significant interest in the su¬ 
perconducting behavior of hydrogen-rich solids at high 
pressures H7HH]. These compounds display a fascinat¬ 
ing variety of crystal structures, as some contain molec¬ 
ular hydrogen, in others the H-H bond is broken, and 
hydrogen forms direct bonds with the other atoms. The 
predicted T c ’s close to ambient pressures are just a few 
K, while for pressures above 150 GPa they range from 64 
to 235 K. 

Among these compounds, H 3 S has been predicted to 
exhibit one of the highest T c ’s. Reference [ 8 ] underlined 


that the mechanism boosting its T c is essentially the same 
one that occurs in magnesium diboride (MgB 2 ), which 
broke the record of critical temperatures for conventional 
superconductors in 2001 j25lf27] . In MgB 2 , the presence 
of electronic bands with a strong covalent character at 
the Fermi level leads to e-p matrix elements much larger 
than those found in other conventional superconductors. 
The experience with this material has shown that, be¬ 
sides the well-known strategy of carrier doping, a very 
promising route to increase T c is to act on the bonding 
properties of the materials, directly influencing the e-p 
matrix elements. Within the range of MgB 2 -like mate¬ 
rials, the highest critical temperature was predicted for 
hole-doped LiBC, whose stiff ionic bonds are even more 
favorable for superconductivity [28 . 

H 3 S is very similar to MgB 2 as it also exhibits strong 
covalent bonds giving rise to large electron-phonon cou¬ 
plings mu and strong anharmonic effects [6l l29ti3T] . and 
although its critical temperature is already spectacularly 
high, there is no reason to expect that it cannot be in¬ 
creased any further. 

The aim of this Rapid Communication is to explore 
whether the superconducting properties of H 3 S can be 
further enhanced by increasing (or decreasing) the cova¬ 
lent character of the sulfur-hydrogen bond. We perform 
a systematic study of the electron-phonon properties of 
several H 3 A compounds in the bee high-pressure struc¬ 
ture of H 3 S, using linear-response calculations [32]. The 
X atom is either one of the first four atoms of the chalco¬ 
gen group (O, S, Se, Te), or a hypothetical alchemical 
mixture of neighboring chalcogens (i.e., O and S, S and 
Se, and so on) [33] [34] . This approximation allows us to 
study the influence of covalency and atomic size on the 
bonding, dynamical stability, and electron-phonon cou¬ 
pling of the H 3 X structure in a controlled way, without 
introducing spurious effects related to doping. We will 
also introduce quantities to measure the bonding char¬ 
acteristics of dense hydrogen-rich compounds that will 
be useful in the search for new high-T c materials with 
high-throughput methods [35]. In order to minimize the 
number of external parameters, we work at a fixed pres- 
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FIG. 1. (Color online) (a) Crystal structure and isocharge 
surface (p = 0.01) of IGSe in the conventional unit cell. The 
Se atoms are indicated by green spheres, and the hydrogen 
atoms by blue spheres. The red surfaces are isocharge surfaces 
(see text), (b) Charge plot in the gray-shaded plane of (a) for 
tCSe and (c) for H3O, both with substracted backgrounds. 


FIG. 2. (Color online) (a) Electronic band structure of H 3 Se 
(solid black) and H 3 0 (dashed red) along the path r(0,0,0)- 
(0,1,0)-iV(0. 5,0.5, 0)-P(0. 5,0.5,0.5)-r(0,0,0)-TV(0.5,0.5,0) 
in units of 27r/a. (b) Partial DOS of H 3 Se with Se (one H) 
contributions in black (green), (c) Partial DOS of H3O with 
O (one H) contributions in red (green). 


sure [36] . 

Figure [I] shows the Im3m bee structure of the H 3 X 
compounds at high pressures: The X atoms (green) oc¬ 
cupy the 2 a Wyckoff position and form a bcc lattice, 
while the hydrogen atoms (H, blue) occupy the 65 Wyck¬ 
off positions in between. Each H forms two short and four 
long bonds with its X neighbors in the (100), (110) and 
(110) directions. When the electronegativity of the X 
atom is slightly larger than that of H ( 2 . 1 ), as in S (2.4) 
and Se (2.5), the charge accumulates around the center of 
the short X-K-X bonds, giving them a strongly covalent 
character. The red, H-centered ellipsoids in Fig. HW de- 
limit regions in which the charge density of H 3 Se is higher 
than that of hypothetical Se and H 3 solids in the same 
crystal structure. A contour plot of the same quantity 
in the basal plane is plotted in Fig. [ljb). As the X-H-X 
bonds acquire an increasing ionic character by increasing 
the electronegativity of X, the maximum of the charge 
density distribution moves from the vicinity of hydrogen 
to that of the X atom; this is seen in the contour plot 
for H 3 0, shown in Fig. [ljc). 

We will now discuss in detail how and why the elec¬ 
tronic and vibrational properties of these two compounds 
differ. We start from Fig. [2|a), where we show the elec¬ 
tronic structure of the two compounds. The band struc¬ 
ture of H 3 Se (solid black) is in nice agreement with pre¬ 
vious calculations and, as discussed by several authors, 
very similar to the one of H 3 S 0018 ]. Se- and H-derived 
electronic states are strongly hybridized, especially in the 
vicinity of the Fermi energy. This leads to several avoided 
crossings and eventually to a strong van-Hove singularity 
slightly below Ep [ 57 ] . 

At the same pressure, H 3 0 would have a lattice con¬ 
stant ^ 25% smaller than H 3 Se; its band structure is 
plotted as dashed red lines in Fig. [2^a). The remark¬ 


able difference in the electronic structure of the two com¬ 
pounds is due to the large electronegativity of O, which 
increases the hybridization between O and H states and 
leads to even larger avoided crossings. As a consequence, 
the van Hove singularity close to Ep splits in two peaks, 
which move apart as well. 

Figure [5] shows the corresponding vibrational proper¬ 
ties: Phonon dispersions are plotted in Figs. [3ja) and 
§c), while Figs, [^b) and |3jd) show the partial phonon 
density of states (DOS) and the e-p spectral function 
a 2 F(uj ), 

O' F(ij) = — ^ ^ l^kjk+qjZ/l ^( e k)^( e k+q)^(^ — <^q,zy) 

Ef kq,^ 

( 2 ) 

where Np F is the DOS at the Fermi level, is 

the phonon frequency of mode v at wavevector q, and 
|#k,k+q,z/| is the electron-phonon matrix element [38], be¬ 
tween two electronic states with momenta k and k + q at 
the Fermi level. The dotted line shows the u;-dependent 
e-p coupling parameter A(c 0 ) = 2 dVt a ; A(oo) 
gives the total e-p coupling parameter already discussed 
in Eq. ([!]). 

In H 3 Se [Fig.[3ja)] the acoustic modes extend up to an 
energy of about 50 meV. The first six optical branches 
are separated by a gap of around 25 meV from the acous¬ 
tic ones and have a bandwidth of ^100 meV. Separated 
again by an energy gap of 25 meV are the last three 
optical branches. These three modes have a mainly 
bond-stretching character, according to the definition of 
Ref. [bj. In turn, the middle six phonon branches have 
a mostly bond-bending oscillatory character. The acous¬ 
tic modes contribute ~ 1/3 to the total value of A. The 
a 2 F(uj) of the bond-stretching modes is distinctly differ¬ 
ent from the shape of the phonon DOS and adds about 
1/4 to the total A. For H 3 Se, we obtain uj\ n = 84 meV 
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FIG. 3 . (Color online) Phononic band structure of (a) ITSe 
and (c) H3O0.5S0.5, where the size of the circles indicates 
the amount of bond-stretching character. a 2 F(cj) (black), 
integrated A (black dotted), and phonon DOS (Se character 
red, H character green) of (b) SeH 3 and (d) H 3 0 . 


and A = 1.35, leading to a T c = 100 K for fi* = 0.1. 

The superconducting properties improve for the 
smaller -Z chalcogens, which have smaller atomic masses, 
decreasing ionic radii, and increasing electronegativities. 
In agreement with previous calculations, we find that 
the T c of H 3 S is distinctly larger than that of Fl^Se (for 
H 3 S we get uj\ n = 97 meV, A = 1.83, and T c = 151 K for 
p* = 0.1) |39|. In addition, we find that a further im¬ 
provement of the superconducting properties occurs in 
the oxygen-rich part of our alchemical phase diagram, 
where the X-H-X bond becomes more and more ionic, 
before the dynamical stability of the lattice breaks down 
eventually for oxygen contents larger than 50%. Indeed, 
H 3 O should crystallize in the Im3m structure only above 
11000 GPa HUE]. 

Figure [3|c) shows the electron-phonon properties of 
the last stable compound, i.e., H3O0.5S0.5. The acoustic 
branches extend up to ^70 meV and contribute again 
about 1/3 to the total A = 1.95. The optical modes are 
located in the energy window from 80 to 245 meV. Bond¬ 
stretching modes are mainly responsible for the peaks in 
a 2 F(uj) centered at 110, 150, and 175 meV. Integrating 
a 2 F{uj) for H3O0.5S0.5, we obtain cji n = 100 meV and 


T c = 164 K for p* = 0.1. 

We have so far discussed only specific choices of X. To 
comprehend how individual differences in bonding prop¬ 
erties, dynamical stability, and electron-phonon coupling 
translate into actual trends in an alchemical phase di¬ 
agram, we present in Fig. [4] an overview of representa¬ 
tive quantities as a function of a continuous variable /3 
that ranges from zero (oxygen) to three (tellurium); non¬ 
integer values represent alchemical pseudoatoms. 

Figure Qa) shows in black the optimized lattice con¬ 
stant, which increases smoothly with f3 from a = 2.62 to 
3.33 A. In red, we plot a compact measure of the charge 
localization along the short bond direction, i.e., the frac¬ 
tion px/ptot of total charge which is localized around 
the X atom [42] . This quantity increases as (3 —)> 0, as 
the maximum of the charge distribution shifts smoothly 
from the H to the X atom. In Fig. [4^b) we show the fre¬ 
quency evolution of the triply degenerate optical phonon 
branches at the Y point - see Fig. [3] All branches have a 
(partially) H bond-stretching and bond-bending charac¬ 
ter and evolve smoothly with /3. While the green and red 
branches have T\ u symmetry, the black one is of T^ u sym¬ 
metry. For high O and Te contents, the lowest branch 
displays imaginary frequencies, indicating a dynamical 
instability of the structure mm- In O-rich compounds, 
another instability occurs earlier near the P point. The 
gray-shaded areas in the figure indicate parts of the (3 
phase diagram where the H 3 X structure is dynamically 
unstable. 

In Fig. [4^c), we plot in red the value of the electronic 
DOS at the Fermi level Ne f ] it is essentially constant 
for S—?>Se—)>Te (/3 > 1 ), and decreases sharply for (3 -> 0, 
due to the increased O-H hybridization - see Fig. [ 2 ] On 
the other hand, A is almost monotonously increasing for 
decreasing /?; uj\ n [Fig. ^e)] exhibits a similar behavior, 
i.e., it remains nearly constant between O and S, and 
bends downwards for Se and Te. A large value of uo\ n for 
small /? is not in contradiction with the incipient lattice 
instability, since, as mentioned earlier, the phonon modes 
which undergo softening carry only a small fraction of the 
total e-p coupling [26, 27j|44]. The combined effect of uj\ n 
and A thus results in a net monotonous increase of T c , 
from 83 K on the Te-rich side to 164 K on the O-rich side, 
as shown as the red dashed line in Fig. §e) ESI 

Figure [4](c) also shows that the behavior of A does 
not correlate with Ne f , especially for /3 0. Instead, 

the average e-p coupling matrix element X/Ne f , plot¬ 
ted in Fig. [IJd), increases by ^ 50% when going from 
F^Se to H 3 O. To rule out that this increase is not just 
a lattice-constant (pressure) effect, we additionally per¬ 
formed e-p calculations for H 3 S at the same lattice con¬ 
stant as H 3 O 0 . 5 S 0 . 5 , i.e., at a pressure of 265 GPa: The 
results, shown as isolated triangles and stars in Figs.^c) 
and[4][d) are in line with those of the other chalcogen hy¬ 
drides (A = 1.70, co>i n = 103 meV, and T c = 152 K). The 
behavior of \/Ne f strongly resembles that of px/ptot , 
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FIG. 4. (Color online) (a) Lattice constant a (solid black) 
and partial charge px / ptot (dashed red). (b) Optical phonon 
frequencies at T. (c) A (solid black) and Ne f (dashed red), 
(d) \/Ne f (solid black) and rescaled partial charge px / Ptot 
(dashed red), (e) T c (solid black) and cj\ n (dashed red). The 
isolated symbols in (c)-(e) indicate results for H 3 S calculated 
at the same lattice constant as H 3 Oo. 5 So .5 (see text). 


plotted in the same panel as the red lines. This is a 
clear indication that increasing the ionic character of the 
X-H-X bond can be beneficial for superconductivity. Al¬ 
though the increase in T c is not as spectacular as in LiBC 
compared to MgE >2 [28], our findings indicate increased 
ionicity as a clear path that experimentalists could follow 
to optimize T c in high-pressure hydrides. 

In summary, we have studied a possible route to opti¬ 
mize the critical temperature in recently discovered high- 
T c superconducting hydrides using first-principles linear- 
response calculations. We constructed hypothetical al¬ 
chemical atoms, which smoothly interpolate atoms of the 
chalcogen group, and studied their phonon and electron- 
phonon coupling properties. The effect of alchemical dop¬ 
ing is qualitatively very different from charge doping or 
pressure, as it allows one to tune the bonding proper¬ 
ties of a material, crucially affecting the e-p matrix el¬ 
ements. We observed the existence of a “sweet spot” 
for electron-phonon superconductivity in H3O0.5S0.5, for 
which the bond is ionic enough to profit from increased 
electron-phonon coupling, but not so strongly as to lead 
to a structural instability. 

This implies that, although H 3 S is already almost op¬ 
timal for electron-phonon coupling superconductivity, its 
critical temperature could be improved further making 
the X-H-X bond more ionic, i.e., partially replacing sul¬ 
fur with more electronegative atoms, such as oxygen 
or members of the halogen group. Together with the 
compact measure we provided to characterize the bond¬ 
ing situation, our results highlight promising paths in 
the search for new conventional high-T c superconductors 
with high-throughput methods [35] . 

Note added. Recently, we became aware of a similar 
investigation [50] . 
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